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1. Introduction 

The Maxwell- Vlasov equations give a model for the dynamics of a large collection of 
charged particles. They are used to analyse the motion of beams of particles in a particle 
accelerator, in order to address problems such as the effects of coherent synchrotron 
radiation and space charge. They are also used for low energy particle dynamics as 
in Klystrons and magnetrons. With several species of particles, the Maxwell- Vlasov 
equations model the dynamics of plasmas both man made, as in laser plasma wakcfield 
acceleration and fusion reactors, and naturally occurring as in the solar winds and the 
ionosphere. 

In this article we formulate the distributional form of the Maxwell- Vlasov equations, 
write down the general solution to the Vlasov equation and identify sufficient conditions 
such that these solutions are valid sources for Maxwell's equations. We then give some 
example solutions. This work is relevant in several areas of active research. 

• It unifies the Maxwell- Vlasov equations and the cold fluid, multicurrent and water 
bag model of charge together with the Klimontovich distribution. 
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• By comparing the statistical limit of the Klimontovich distribution with the regular 
solutions of the Maxwell-Vlasov equations, it may provide a way of finding a 
dispersion term for the Boltzmann equation. 

• In many scenarios, such as micro-bunching and emitance reduction in accelerators 
and charge moving on surfaces of constant magnetic flux in a Tokamak, the 
resulting charge distribution may be better modelled by the use of lower dimensional 
distributions. 

• It enables the description of the ultra-relativistic expansion for the Maxwell-Vlasov 
and the Klimontovich distribution. Since a single system applies to all the solutions, 
finding the expansion of this system will enable one to write down the expansion 
of all the solutions above. 

• It enables the use of both the retarded Greens potential and the Lienard-Weichart 
potential to find the electromagnetic field due to an arbitrary distributional source. 

• It enables alternative methods of finding numerical results. The usual method is to 
use a collection of world-line distributions based on the motion of individual point 
particles. Alternative methods could use ribbon or higher dimensional distributions. 
This may avoid some of the regularisation problems associated with point charges. 



1.1. Notation 



Sections [2] and [3] deal with general properties of distributions and the transport equations 
(also known as the Vlasov, Liouville and coUisionless Boltzmann equation). For generic 
objects in these sections we use the following symbols: 



Manifolds: M,N,P,Q. 

Subsets: U,V C M. 

Boundary: dM 

Tangent bundle: TM. 

Bundle of p— forms: A^M. 

Bundle of forms: AM. 

Generic bundle: it : E ^ M. 

Smooth sections of a bundle: TE. 

Space of test forms: TqAM. 

Space of continuous (not differentiable) 

forms: FctsAM. 
Space of piecewise continuous 

forms: FpcAM. 
Space of test forms: FqAM. 
Space of distributions: FdAM. 
Set of submanifold distributions: FsAM. 
Vectors and vector fields: u, v, w. 
Forms and form fields: a, f3. 



Evaluation of a field at a point: a\x, u\ 
Test forms: 0, ip. 
Distributional forms: \E', $. 
Regular distributions: D{a). 
Degree of a form or distribution: 

deg(a), deg(^'). 
Sign involution a"^ = (— 1)'-^^S(")q,. 
Smooth map between manifolds: 

a, b, c, a: ^ M. 
Image set of the map a(A^) C M 
Preimage set of the subset U C M: 

a-\U) C N 
Embedding: a: N M. 
Composition of maps: a o b. 
Pushforward for distributions: a<j, b<^. 
Submanifold distribution: a^a. 
Fullback for distribution: a''. 
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Box: s:S^M. Curve: 7: M M. 

Expectation of a distribution with Lift of a curve: 7 = 7^(9^): M "—>■ TM. 

respect to box s: S M: [s''(\l/)]. Initial hypersurface: a: S ^ M. 

In sections [T] and |4] we look at the Maxwell-Vlasov equation on spacetime. We use 
the following symbols: 

Spacetime: A4. 
Spacetime metric: g. 
Upper unit hyperboloid: £ 
Metric dual of a 1-form: 5, 
Electromagnetic 2-form: F 
Liouville vector field: W. 



Charge distributional 6-form: B. 
Electromagnetic source distributional 3-form: J'. 
Worldline: C:R^M. 
Lift of the worldline: C = C^(9^): M ^ £. 
Double hft of the worldline: C = ^(9,): M ^ TS. 



1.2. Maxwell-Vlasov equations 

Let {M.,g) be spacetime and let 

Ti:£^M where £ = [u eTM\g{:u,u) = -I, > (1) 

be the unit upper hyperboloid bundle over spacetime. Let the coordinates 
(x", X"*^, x^, x^, y^, y^), for £ with the embedding £ ^ TM. be given by x" = x", 
X* = and let x° = ?/° = ?/°(x",?/*) be the solution to gabU'^y'' = —1, 2/° > 0, and 
yo = QoaV"'- Here the indices a, 6, c = 0, 1, 2, 3 and i, j, k = 1,2, 3. We choose physical 
units of time and length so that the speed of hght c = 1 and the permittivity of free 
space eo = 1. 

Given the electromagnetic 2-form field F G FA^A^, the standard way of expressing 
the Maxwell-Vlasov equations is to prescribe the Vlasov vector field W = W{F) G TT£, 
which depends on the electromagnetic field as 

w = y''^a- y''f^, + -y"" g'"^, (2) 

for the species of particle with mass m and charge q. For most of this article we deal 
with a single species and therefore choose physical units of mass and charge so that 
m = 1 and q = 1- 

This vector field is chosen so that is it horizontal, i.e given u & £ then vr^(W^|u) = u 
and that if 7: M — > £^ is an integral curve of W, that is 7(r) = 7*(9t-) = then 
7 = (7 = Ci,{dr) where C = 7ro 7: M^^Vl satisfies the Lorentz force equation 

VcC = (3) 



with g{C,C) = -1, see [H |2] and lemma 29, section 3.1 Here~: A^M TM is the 



metric dual given by l3{a) = g a). A function / G rA^£ is called the one particle 
probability functionj|] The Vlasov equation, which is also know as the Liouville 

I The use of the word distribution is avoided in this context. This word will be reserved for distribution 
in the sense of Schwartz or De Rham currents. 
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equation and the coUisionless Bohzmann equation, is now given by 

W{f) = (4) 
Maxwell's equations for the electromagnetic field F are given by 

dF = and di<F = -J (5) 
where the source J G FA^TVI is given by 

J=( [ — V\d^\dy''%/9x^ * 1 (6) 

where dy^"^^ = dy^ A rfy^ A dy^ etc. There is a natural non-vanishing 7-form on S given 

by 

yo 

Since the Lie derivative L^yQ = then 

dtwim = Lwim = Lwim + fLw^ = Lwim = wim 

Thus Q is equivalent to diwif^) = 0. 

We introduce the charge 6-form given by ^ G TA^S 

9 = lwim (8) 

Thus we can rewrite Q as 

de = and iwO = (9) 

Observe that given 6 and Q it is easy to construct the probability function / since the 
coefficient of 9 with respect to dy^"^^ A dx^^^ is given by f det{g). 

Comparing ^ and ([s]) we see that the integration for the source G TA^Ai, is in 
fact a integration along a fibre [3] given by 



7r>Ae= / 0AJ (10) 

for all test 1-forms G FqA^AI, that is the space of smooth 1-forms with compact 
support on A4. 

The 6-form 6 has a natural interpretation in terms of probabilities or charge. Given 
a "box", that is a compact 6-dimensional region S G S which is transverse to W 
then the probability that the particle passes though S (or the total charge passing 
though S) is given by the integral jg 9. An example of such a box is the set S = 
{(x", = , < < , yl < y^ < Vu} for some constants {x^^^x\,x\,y\,y\^}. 

Assuming the correct orientation for S it is usual to demand j^O > 0. An initial 
hypersurface S C £^ is any 6 dimensional hypersurface transverse to W. If 6 is to be 
interpreted as probability then we further demand that Jj,6 = 1. However if 6 is the 
charge distribution then = Qtotai is the total charge. It is easy to show that Qtotai 
is a conserved quantity. In order to evaluate a distribution on a box or a hypersurface 



we require the definition of the puUback of a distribution given in section 2.3 
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In order to write the Maxwell-Vlasov system in distributional language we replace 
the 6-form 6 with the charge distributional 6-form 

QeTj^A^S, (11) 

This is defined as an element in the dual to the space of test 1-forms TqA^S. Given a 
test 1-form G TqA^S then the action of on is the real number written 0[0]. 



Inspired by (10) we define the pushforward of a distribution vr^ [4j via 

7r,(e)[0] = 0[7r^(0)] (12) 
and we have the source for Maxwell's equations as 

J e r^A^M ; J = 7T,Q (13) 



The details of this are given in section 2.1 



Since J" is a distribution, then in general, the solution to Maxwell's equations ^ 
will also be a distribution. However in order for F to drive the Liouville equation ( 2p ) 



we require that F is continuous, though not necessarily differentiable. Since F is not 
differentiable, we must convert it into a distribution before inserting it into Maxwell's 
equations. Every continuous g-form a G FctsA'^M gives rise to a regular distribution q- 
form D{a) G TdAW. Thus we can define D{F) G TdA^A^ given by D{F)[(p] = 0AF 
for any (p G ToA^Ai. Maxwell's equations become 

d{DF) = and d{D{^F)) = -J (14) 

As is noted in section |2] we can define the exterior derivative and the internal contraction 
of a distribution, thus the Liouville equation ^ becomes the equivalent for namely 

de = (15) 

and 

iw'S) = (16) 
where the Liouville operator W is given by ([2]). These are also called the transport 



equations of G for W. Theorem [27| implies that the source is closed, dj' = 0. 

We limit ourselves to "submanifold distributions". These are written G = a^a 
where a: £^ is an an embedding with possibly having a boundary and a G FA^A^ 
is a smooth p— form where p = dim A^— 1. The action of this distribution on a test form 
(p G TqA^S is given by 

a^a[(f)] = / a* {(f)) A a (17) 
Jn 

Most of section |2] details the mathematics needed to manipulate these distributions. 

The geometric nature of these distributions means that we can define many useful 
properties such as puUbacks, without having the worry about the convergence and 
other topological properties. However these distributions are sufficiently general that 
physically interesting solutions to the Maxwell-Vlasov system can be written using 
them. In section |4] we show that the cold fiuid, the multicurrent and the water bag 
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models of charge dynamics are all examples of submanifold distributional solutions to the 
Maxwell-Vlasov equations. Another example of a submanifold distribution, discussed 
in section 4.1, is the worldline for a single charged particle given by = C,;(l) so that 
G[0] = JjgC*(0) for all e TqA^S. In Minkowski space, this is equivalent to the usual 
expression for a charged particle, written in terms of a Dirac 5-function probability 
function as 

3 3 

/ = ( n S{x' - C\r))) ( n S{y^ - C'ir))) 

i=l j=l 
where t = C^{t) so that r = {C^)^^(t). Thus we have 

c,ii) = iwim (18) 

To see this recall that C : M ^ £^ is an integral curve of W so that W^|c'(t) = C(r) 
where C = C^,{dr), thus given a test 1-form (f) G TqA^S we have 

(7*01, =idr{C*^)\rdT =C''{ic^^g^y)(P)dT = C^ic^^^cfydj =C* {iw\^^^^4>)dT = {iw4>)\c[r)dT 

Hence 

Q(1)[0]= / C\<P)= [ M)\cir)dr 



and 



tw{m[<i>] = imi^w 

3 



[^w(t>) 5 



[iw4>)\c{r)l^;^dt = / {iw(t))\c^^)—dt = / {iw(t>)\c{r)dr 



Equation (18) is true for general spacetimes, assuming the correct definition of the 
^-function. One of the advantages of the pushforward approach to distributions is 
that the definitions do not explicitly contain the coordinates system. They also do not 
depend on the metric and therefore can be applied to the upper unit hyperboloid bundle 
£ without having to choose a metric for this manifold. 

The Klimontovich distribution is given by the sum of worldline distributions one 
for each particle [5] . We show that the worldline of a point charge and the Klimontovich 
distribution satisfy the Liouville equation and have the correct form for the source of 
Maxwell's equations. However the solutions to Maxwell's equations with the worldline 
as a source, which in free space are given by the Lienard-Weichart potentials, are not 
continuous and therefore we cannot insert them into the Liouville equation without 
some form of regularisation, leading, for example, to the Lorentz-Dirac equation. The 
same is true for the Klimontovich distribution, where the usual regularisation is simply 
to ignore the field due to a particular particle when looking at the dynamics of that 
particle. 

In section |3] we look at the submanifold distributional solutions to the transport 



equations for a general vector field v E TTM on a general manifold M. Theorem 21 
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Figure 1. A submanifold solution to the transport equation. Observe that the 
submanifold is tangential to the vector field and that the boundary of submanifold 
is also tangential to the vector field. 



gives the necessary and sufficient conditions for a submanifold solution. A sketch of 
such a solution is given in figure [TJ We look at initial hypersurfaces and boxes and 



show the conservation of charge, theorem 22, as well as that, given a distribution on an 
initial hypersurface, there exists a unique solution to the transport equations, theorem 
24 This is achieved by extending the initial value along the integral curves of v. 

Given a solution B to the Liouville equations on S we can ask if J' = tt^Q is a 
valid source for Maxwell's equations. The problem is that, in general, tt^Q may not be 
a distribution. This is because vr is not a proper map. A map is proper if the preimage 



of any compact set is compact. This is required so that 7r*((/)) in (12) has compact 



support. In section 3.1 , lemma 28, we show that an appropriately bounded submanifold 
distribution is a valid source. We say that is bounded with respect to vr. This 
is the case in the examples given in section |4j We then pose the following question 



(70): If a distribution is initially bounded, does it remain bounded? Surprisingly the 



answer to this question is no and two counter examples are given, one using a drifting 
source in Galilean spacetime, the other, on Minkowski spacetime, where the initial data 
is specified on a lightlike surface. However if, as we generally assume, is a globally 
hyperbolic spacetime, and the initial data is given on a Cauchy surface, then theorem 30 
shows that the solution will remain bounded. However, as in the case of the worldline for 
a point charge, the source may be valid distribution but the resulting electromagnetic 
2-form need not be continuous. 



It is easy to extend (13 16) to model several species of particle with masses 
and charges qa- Let Bq G FdA^^^ be the corresponding one particle 6- form distribution, 
and Wa G TT Ai be the Vlasov vector field ^ with appropriate m and q. The Vlasov 
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equations are given by 

dQa = and iw^Qa = (19) 



and Maxwell's equations are given by (14) where the source 



:^ = $^vr,(0J (20) 



2. Distributions 

Let M be an arbitrary manifold of dimension m. Let FqA^'M be the space of test 
]>forms. 



FnA^M 



|0 G FA^M has compact support | (21) 

Let the space of p-forms distributions be the vector space dual of the space of test 
functions. 

FdA^M = (FoA'"-PM)* (22) 

We use square bracket notation. 

FdA^-^M X FoA^M ^ M , (^, 0) ^ $[0] G M (23) 

Observe that the space of piecewise continuous p— forms is a subspace of the space of 
p— form distributions. These are called regular distributions. 

D: FpcA^M FdA^M ; D{a)[^] = [ (pAa for (peToA'^-PM (24) 

The exterior derivative of a p— form distribution is given by 

ci: FdA^M ^ FdAp+^M ; rf^[0] = -^[^0"] for G FgA'^-P+^M (25) 

The internal contraction of a vector field v G TTM with a p— form distribution 

G FdA^M is given by 

TTM X FdA^M ^ FdAp-^M , {v, ^) i^^ (26) 
where i^^[0] = -^[^^,0''] for all G FqA'^-p+^M 
Given a regular distribution -D(/3) with /? G FpcAM then 

z,(D/3) = D{i,f3) (27) 

Given a regular distribution D{f3) with /9 G F^tsAM and (3 has a continuous derivative 
and where M has no boundary then 

d{D(3) = D{d(3) (28) 

We define the support of a distribution «S(^) C M as 

5(^) = p||t/cM f/ is closed and if G FqAM, 5(0) C M\f/ then ^[0] = o} (29) 

where «S(0) is the closed support of 0. That is «S(\E') is the smallest closed subset of M 
such that for every test form with «S(\E') fl «S(0) = then \E'[0] = 0. 
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Lemma 1 For regular distributions S{D{a)) = S{a). 

proof: Clearly if 5(0) H S{a) = then D{a)[(p] = / A a = 0. Thus S{D{a)) C 5(a). 
If S{D{a)) 7^ 5(a;) then there exists an open set U C S{a)\S{D{a)). It is easy to 
construct a test form with 5(0) C U such that J A a 7^ 0. So S{Da) H U ^ 
leading to the contradiction S{Da) fl (5(a)\5(Da)) 7^ 0. Hence result. ■ 

Lemma 2 For ^ G TdAM and G TTM 

5(d^) C 5(^) and 5(z„^) C 5(^) (30) 

j)roo/.- Given G TqAM such that 5(0) C M\5(^) then 5(c?0) C M\5(^) and 
5(i^0) C M\5(^). Thus = -*[d0^] = d^[0] and = -^[z^0''] = ^,,^[0]. Hence 



M\5($) C M\5(d$) and M\5($) C M\5(i„$) and (30) follows. 



^. i. Pushforward of distributions 

The pushforward |^ of a distribution with respect to a smooth map a: N ^ M where 
dim M = m and dim = n is given by 

a,: TdA^AT ^ TdA^-^+pM ; a,(^)[0] = ^[a*(0)] for G ToA^-^M (31) 

The pushforward does not preserve the degree of a distribution. Instead we have 

deg(a^(^)) = deg(^) + dim M - dim (32) 



We note that the pushforward given in definition ( |31[ ) is not always defined. The 
problem is that in general there is no guarantee that a*(0) has compact support. This 
can lead to the pushforward being infinite as in the following example. Let a: M — >■ {x} 
and let = D{dt) G rDA^{2;}. For G roA°{x} we have 

a^(^)[0] = ^[a*(0)] = / 0(x)dt = 0(x) / dt 



which is undefined if 0(x) 7^ 0. It can also lead to problems with boundaries. 

By contrast, we can define the pushforward of a distribution when the map 
a: A^ — > M is proper, that is, when the preimage of every compact set is compact. 



However the bundle map n: S ^ M. is not proper but from (12) we see we wish to 
define the pushforward of a distribution with respect to n. This is not always possible, 
but it is for certain distribution which we say are bounded with respect to vr. 

We say the distribution \E' G r^AA^ is bounded with respect to a: A^ — ^ M if for 
any compact set ?7 C M then 5(^) fl a^^{U) G N is compact. In this case we define 

a^(^)[0] = ^[/ia*(0)] for G TqAM (33) 

where h G FqA^A^ is any test bump function such that 

5(*) n5(l - /i) n5(a» = (34) 

Lemma 3 The pushforward of a distribution bounded with respect to a is well defined. 
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proof: Since is compact and «S(a*0) C a~^(«S(0)) then S{a*(f)) n«S(\I/) is compact. 
There exists a compact set Vi and V2 such that «S(a*0) fl C interior(l^i) and 

V2 C interior(Vi) and a bump function /i such that h{x) = 1 for x G Vi and = 
for X ^ V2. 

To see that a^(\l/)[0] is well defined, consider two functions h and h, then since 
h-h={l-h)-{l-h)we have S{h - h) C 5(1 - /i) U 5(1 - h). Thus 

5(^) n5(a» n5(/i-/i) 

C (5(^) n 5(a» n 5(1 - /i)) U (^5(^) n 5(a*0) n 5(1 - /i)) = 

Thus ^[(/i-/i)a*(0)] = 0. ■ 
To guarantee that the pushforward is well defined we limit ourselves to three cases: 

• When a: — > M is proper. 

• When \1/ is bounded with respect to a. 

• When 7c: E —>■ M is a. fibre bundle and a G TAE has the property that integration 
along the fibres tt~^{x), given by 7r*0 A {n^{Da)) = /^^ A a for all is defined 
[3]. This last case is simply for regular distributions which are not bounded with 
respect to vr but such that integral exists. We will not consider this standard type 
of pushforward further in the article. 

Lemma 4 Given manifolds N, M, P and smooth proper maps b: P N and a: N ^ M 
then the composition of the push forwards is the push forward of the composition: 

(a o b)^ = a^ o (35) 

proof: The composition of two proper maps is proper. Let \1/ G FqAM and (f) G FqAP 
then 

a,(b,(^))[0] = b,(^)[a>] = ^[b*a*0] = ^[(ao b)>] = ((ao b),^)[0] 



Observe that if a or b is not proper then care must be taken. See lemma 28 



Lemma 5 If a: N ^ M is proper then the pushforward of distributions commutes with 
the exterior derivative. 

do a^ = a^ o d (36) 

proof: 

d{a^^)[(f)] = -a,^[(i0^] = -*[a*(t^0'')] = -^[rfa*(0'')] = -^[rf(a>^)] = rf^[a*(0)] 
= a,{d<fM 

m 

Lemma 6 Let a: N ^ M be smooth and let \E' be bounded with respect to a, then 

(i(a^*) = a,(d*) (37) 
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proof: From (30), d"^ is bounded with respect to a. Given h e ToA'^N such that (34) 



holds, then from (|30]) n 5(1 - h) n S{a*(j)) = and since S{dh) C 5(1 - h) we 

have 5(^) n S{dh) n 5(a*0) = so "^[dh A a>^] = 0. Thus 

d{a,^)[(j)] = -a.^ic/^"] = -^[/ia*(rf0^)] = -^[/i rfa*(0^)] 

= -*[d(/iaV)] + ^[dh A aV"] = -^[rf(a>^)] = (rf^)[a*(0)] = a,(rf^)[0] 

■ 

Given a: N ^ M and a vector field we say that v G TTM is tangential to a if there 
exists u G FTA^ such that for each x E N, a-^{u\x) = f |a(x)- We write a^('u) = v\a{N) if 
a^(M|a;) = f |a(x) for all x e N, 

Lemma 7 Let a: N M be proper with v G TTM tangential to a and let u G VTN 
satisfy a*(M) = f |a(Af) then 

iy o a,; = a^ o i„ (38) 
proo/; Given \E' G TdAM, G TqAM, a: N ^ M and x G A^, let ?/ = a(x) then 

(i„a*0)|^ = ^„|,a*.(0|y) = a^(ia.H.)(0W) = = ^liM)\y) = {a'' {iv(j)))\x 

where a*: Aj,A^ — ^ A^-M. Thus i^a^^ = a^ii^cj)). Hence 

2„a,(^)[0] = -a,{^>)[i,<P''] = -■^[a\i,(t>'')] = -^[^^a^"] = z„^[a>] = a,(z,^)[0] 

■ 

Lemma 8 Let a: N M be proper and G FdAA^ then 

5(a,^) C a(<S^) (39) 

proof: Given any x G a^^(M\a(«S ^f)) then a{x) G M\a(«S \E') and so a(x) ^ a(«S \1') so 
X ^ 5(^) hence x G A^\5(^). This gives a-^(M\a(5^)) C A^\5(^) 
Let (f) G TqAM with 5(0) C M\a(5 ^). Then 

5(a» = a"^(50) C a~^(M\a(5^)) C A^\5(^) 

Thus a,^[0] = ^[a>] = 0. 

Hence we have shown 5(0) C M\a(5\E') implies a^\E'[0] = 0. Hence M\a(5\l') C 



M\5(a,^) and hence (39). 



To see that in general 5(a^\E') 7^ a(5\l') consider the following counter example. 
Let M = M and = M X {-1,1} with a: A^ M given by a{x,i) = x. Let 
\E' = D{a) G FdAA^ where a G FA^A^ is given by = 1 and q;|(^^_i) = — 1 then 

given G FqA-'^M we have 



a?a[0] = / a*(0)Aa= / 0U A + / 0UAa|(^_i)= / 0- / = 



Hence 5(a,^) = ^ M = a(5 *). 
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2.2. Submanifold Distributions 

A special kind of distribution is the submanifold distribution. Before we introduce 
submanifold distributions we observe some basic properties of closed embedding, i.e. 
embedding a: N M where a(A^) C M is closed. First observe that a is proper. 

Lemma 9 Given an embedding a: N ^ M such that a{N) G M is closed it follows that 
a is proper. 

proof: Given a compact U G M then a(A^) fl f/ is closed and therefore compact. The 
restriction a|a-i(^): a^^(t/) Si{N) fl [/ is a homeomorphism so a^^{U) is compact. 
Thus a is proper. ■ 
The next three lemmas use the embedding nature of a to relate a^{Da) and a. 

Lemma 10 Let a: N ^ M be an embedding. Given an open set U G M and G TAN 
such that «S(0) C a^^{U) then there exists ip G TAM such that a^ip = and Slip) G U . 

proof: Consider x G a{N). Since a is an embedding there exists an open coordinate 
patch V G M, (x^, . . . about x and a coordinate patch a~^V G N, {y^, . . . ,y"-) 
such that a{y^, . . . ,?/") = (x^, . . . ^x'^^O, ... ,0). 

Assume for the moment that «S(0) C a~^{V) and (f) = (pidy^ where / C {1, . . . n} 
and dy^ = dy^^ A ■ ■ ■ A dy^\'\ refers to multi-index notation, then let 

^= J2 (pidx^h{x''^\...,x"') 

/C{l,...n} 

where h{x^~^^, . . . , x™) is a smooth function with h{0, . . . , 0) = 1 and so that S{ip) G V. 
Thus a^ip = (j). 

Now in general use a partition of unity to partition U into coordinate patches Vi. 
Thus if = ^ 0j with «S(0j) C a~^Vi then ip = "^ipi and a*ip = 0. ■ 

Lemma 11 Let a:N^M be an embedding with a{N) G M closed then 

S{a,{Da)) = a{S{a)) (40) 

proof: Lemma [l] implies S{Da) = S{a) so a{S{Da)) = a(«S(a)). Lemma [s] implies 
S{a,{Da)) G a{S{Da)). Thus S{a,{Da)) G a(5(«)). 

By contrast let x = a{y) G a(«S(a)) so y E S{a) as a is injective. Given 
any neighbourhood U G M of x then a~^{U) G N is a neighbourhood of y. As 
y G S{a) there exists ip G FqAA^ with «S(0) C a~^U such that Jj^ (f) A a ^ 0. 
From lemma 10 there exists a ip E FqAM such that «S('0) C U and a*^p = 0. Thus 
ai;{Da)[ip] = Da[a*'0] = Da[(j)] = Aa 7^ 0. Since this is true for all neighbourhoods 
U about X then x G S{a^{Da)). Thus we have shown x G a{S{a)) =^ x G S{a^{Da)) 
hence (|40]). ■ 

Lemma 12 If a: N M is an embedding, a{N) G M is closed and a G FAA^ then 
a^{Da) = if and only if a = 0. 



proof: Follows trivially from lemma 



11 
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A submanifold distribution is a distribution of tlie form \1/ = a^{Da) where 

• a: M is an embedding. (41) 

• ae TAN with 5(a) = A^. (42) 

• a(A^) C M is closed. (43) 

• a: N a{N) is a diffeomorphism. (44) 

The set of all submanifold distributions over M is written FsAM. 

Since the combination of a pushforward of a regular distribution, for example 
a^{D{a)) is so common, we introduce the notation (the bold subscript <;) to represent 
the pushforward of a regular distribution so that a^a = a^{D{a)). 

Observe that if a: A^ M and a^a G FsAM is a submanifold distribution then 



from lemma Pais proper and from lemma 1 1 



S{a^a) = a{N) (45) 

Given \E' G FsAM, the following lemma establishes the essential uniqueness of the 
embedding and the form on the domain of the embedding. 

Lemma 13 Given two submanifold distributions a^a G FsAM and b^(3 G FgAM with 
a: N ^ M and b: P > M then a^a = b^[3 if and only if there exist a diffeomorphism 
c:N-^P with a = c*P and the following diagram commutes 

N 





M 

proof: Since S{a^a) = S{b^l3) then a(A^) = b{P). Furthermore a: N ^ a{N) and 
b: P ■—>■ b{P) are diffeomorphisms. Thus we can let c = (bj^^p^"^ o a: A^ P, so c is a 
diffeomorphism and a = bo c. 



Given G FqAP, from lemma 10 there exists ip G FqAM such that b*ip = 0. Now 



0A/5= / b*'^ A /3 = b^{f3)['ilj] = a^{a)['ilj] = / a*V^ A a = / (boc)>Aa 

Jp J N JN 



c*b»Aa= / c>Aa= / {c* (f) A a) = / A (c~'*a) 

N JN Jp Jp 

Since this is true for all then j3 = c~^*a. ■ 
The following lemmas relate embeddings with internal contraction. 

Lemma 14 Let a:N^M and a^a G FsAAf. Let x G a(A^) and let v G TTM such 
that v\x ^ Tx{a{N)), that is v is transverse to a at x. Then ^^(a^a) 7^ 0. 

proof: There exists an open neighbourhood U C M of x such that v\y is transverse to 
a^a for all y G a(A^) fl U . Shrinking U if necessary we can assume that U fl a{dN) = 0. 
Further shrinking U we can make it contain a coordinate chart about x adapted to v. 
Thus there exists a t:U — * M such that t{y) = for all y G a(A^) and v{t) = 1. 
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Now for ip G To AM with S{ijj) C U then 

a^a[dtAij]= / a* {dt A ip) A a = / da* (t) A a* (ip) A a 
Jn Jn 



a*{t) A a*{^p) Aa- a*{t) d{a*{ij) A a) = 
dN Jn 

since a*t{z) = t{a{z)) = as a{z) G a{N), and since also a*ip\dN = 0. 

There exists (p ^ TqAM such that C U and a^a[</)] 7^ 0. To see this consider 

the contrary that S{(f>) C U imphed a5a[0] = 0. This would imply that U C M\S{a^a), 
a contradiction. 

Now consider iv{a^a)[(f) A dt] 

iv{a^a)[(j) A dt] = —iv{a^a)[{dt A 0)''] = a^a[iv{dt A 0)] = a5a[0] — a^a[dt A i^cj)] 
= a^a[(j)] 7^ 

Hence 2„(a^a) 7^ 0. 

■ 

Lemma 15 Let a^{a) e TsAM with a: N ^ M . Then 

i^a^a = (46) 
if and only if v is tangential to a and 

iua = (47) 
where u G TTN is the unique vector field satisfying a^(u) = t>|a(Ar)- 
proof: That (47) implies (46) follows trivially from lemma [7[ 



If ([46|) is true and x G a(A^) then from lemma 14 must be tangential to a. Then 



lemma [7| gives a^{iu(y) = 0. Finally (47) follows from lemma 12 ■ 

The set of p— form distributions FdA^M forms a vector space in that it is closed 
under addition and multiplication by a scalar. Also the exterior derivative d: FdA'^'M — »• 
TdA'P'^^M. By contrast the set of p— form submanifold distributions FsA^M does not 
in general form a vector space. For example if a^a G FsA^M and b^f3 G FsA^M 
with a:N^M, b: P ^ M and dim ^ dim P then a^a + h^(3 ^ T^A^M. Also if 
dim A^ = dimP and a(A^) ft b(P) ^ then in general a^a + b^/? ^ FsA^M. 

Likewise in section 2.4| we see that the exterior derivative does not map FsA^M to 
FsA^+^M. This is a due to the possible additional boundary terms. 

2. 3. Pullback of Distributions 

We first define the pullback for a general "pushforward" distribution. This is a 



distribution of the form a^a G FdAM. This is required for the proof of lemma 35 
However in most cases the map a is an embedding, and this simplifies the concept. This 
definition of a pullback of a distribution can be compared to the definition using the 
weak limit [U]. 
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Let a: N —>■ M and b: P M he proper maps and let a G FAA^. Thus a^a G FdAM 
is a distribution and we wish to define the pullback b''{a^a) G FdAP with respect to b. 
Let 

Q = {{y,p)eNxP\a{y) = b{p)} (48) 

be the induced manifold (also know as the pullback manifold), such that the following 
diagram commutes 



Q 



(49) 



where b{y,p) = y and a{y,p) = p. 

We say that a^a G FdAM is transverse to b if 

dim Q = dim + dim P - dim M (50) 
Given a^a G FsAM transverse to b then we can define the pullback of a^a by b 

b\a^a) = a^ilj'a) G FdAP (51) 
Lemma 16 The pullback preserves the degree. 



proof: From (32) we have 



deg(b^(a5a)) = deg(a5(b a)) = deg(b a) + dimP — dimQ 
= deg(Q;) + dimM — dim A^ = deg{a^a) 



Lemma 17 For regular distributions we observe 
a'iDa) = D{a''a) 



(52) 



proof: For regular distributions (49) becomes 



P 



Id 



P 



M 

Id 

■M 



Hence for a G FAM, a'-{Da) = Id^(a*a) = D{a*a). ■ 

Lemma 18 Composition of pullbacks: Let b: P M and c:Q ^ P, where b"" and c"^ 
are defined, then 

(boc)^ = c^ob^ (53) 

proof: Let a^a G FqAM with a: N ^ M and a G FAA^. Define P, Q, a, b, c, a using 
( 48p9 ) so that the following commutes 

b 
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then dimQ — dimQ = dimP — dimP = dimA^ — dimM, so that the pullbacks are 
defined. 

(b o c)'"(a^a) = a^((b o c)*a) = a^(c*(b*a)) = c^(a^(b*a)) = c^(b^(a^a)) 



hence (53). 



It is often the case that both a:N^M and h:P^M are embeddings. Observe 
that in this case the intersection Q = a{N) fl b(P) since 

a(b(Q)) = b(a(g)) = a{N) f] b(P) 

and all these maps are embeddings. 

2.4- Boundaries 

Lemma 19 Let N be a manifold with dimA^ = n and have boundary B = dN with 
l: B ^ N and let a G FA^A^ (a is smooth). Then 

d{Da) = D{da) + L^{i* a) (54) 

proof: Let G FqA^^^'-^A^ then 

d{Da)[^] = -Daldcp"] = - [ d(f)'^ Aa = - [ 6/(0" Aa)+ I (p A da 

Jn Jn Jn 



L*{(j)'^ Aa) + D{da)[(j)] = D{da)[(j)] + (-1)"-^ / i*(0) A L*{a) 

'B Jb 

D{da)[(f)] + {-l)''~PD{L*a)[L' 



We observe that in general the right hand side of (54) is not a submanifold 
distribution, since the domain of D{da) is A^ whereas the domain of L^{L*a) is B and 
dim A^ 7^ dim P. 

Also in general L^{L*a) is not a submanifold distribution. This is because in general 
l: B ^ N is not an embedding. For example, if A^ is a 3-dimensional solid bounded 
cylinder then P consists of two discs and a cylinder x I where / is a closed interval, 
and l: B ^ N is not injective. 

Although this doesn't matter for smooth forms since the set where t is not injective 
has measure zero, one has to be careful when dealing with distributional forms. 

Lemma 20 Let a G FAA^, so that a is smooth then 

d{Da) =0 =^ da = (55) 

proof: Let P = dN with t: B ^ N, deg(a) = p and dim A^ = n. Given G FqAM with 
5(0) C A^\P then 

= d{DaM = D{daM + (-l)"-''i,(i*a)[0] = P(rfa)[0] + {-iy'~fD{i*a)[L*(P] 
= D{da)[(j)] 
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since L^cj) = 0. Thus S{D{da)) C B. From lemma [T| S{da) C B. So da = on N\B, 
but since da is continuous and B contains no open sets, da = 0. ■ 
If a: N ^ M and neither M nor has a boundary then we have the trivial result 



following from ( 28 ) lemma [5j 
d{a^a) = a^{da 



(56) 



3. Distributional solutions to transport equations 

Given a manifold M, with dimM = n and a nowhere zero vector field v G TTM then 
we say that the distribution \E' G rj)A^~^M is a solution to the transport equations if 

rf^ = and = (57) 

We use the name transport equations to refer to any vector field v G TTM. The 
transport equations with respect to the vector field W G TT£ where tt-.S^ A4, 
£ C TM is a bundle (usually the upper unit hyperboloid over spacetime) and W is 
horizontal, that is 7r^(Vr|u) = u for u ^ S, are also known as the Liouville equation, the 



collisionless Boltzmann equation and the Vlasov equation (15 16). 



In general we will only be interested in solutions to the transport equations which 
are submanifold distributions and which, if they possess a boundary then the boundary 
is an embedding. See figure [Tj These are given by 

Theorem 21 Let N be a manifold which is either without boundary or if it possess a 
boundary B then let u B ^ N be a closed embedding and let M be without boundary. Let 
a:N^M and a G rA'^'"'^"^M. Then the submanifold distribution = a^a e TsAM 



satisfies the transport equations (57) if and only if the following three conditions hold 

• da = (58) 

• There exists u G TTN such that a*(M) = v\a(N) and i„a = 0. (59) 

(60) 



If N has a boundary then u is tangential to l. 



proof: First show that the three conditions (|58)|60|) imply that is a solution to the 

then 



transport equations (57). Let w G TTB satisfy l^w 



u\ 



i{B) 



L [Ina] 



6*0 = 



Now since deg(t*a;) = dimA^ 
From lemmas H] and | 



1 = dim(i?) and w 7^ then t*a = 0. 



d^ = da^a = a^((i(Da)) = a^{da) — a^L^{L*a) 
and from lemma [7] 

iy"^ = iya^a = a^{iu{Da)) = a^{D{iua)) = 



In order to show that the transport equations (57) imply (58 60): From (57) we 
have = iy"^ = iya^a. Hence from lemma 15, f is tangential to a and we have (59). 

Let l: B N he the boundary of N. Then = d{a^a) = a^{d{Da)). From lemma 
12 = d{Da)) hence from lemma 20 da = and hence a^L^{L*a) = 0. Since aot: B ^ M 
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is a closed embedding and = a^t^(6*Q;) = (a o L)^{L*a), it follows from lemma 12 that 
L*a = 0. 

Now a satisfies the transport equations with respect to u, i.e. da = and i„Q; = 0. 
If u is transverse to l at any point in B then u is transverse to l in an open subset 
V G B. Furthermore since L*a = 0, iu<y = and u is transverse to l then a\y = for all 
y E V. Using the transport equations for a this implies there exists a an open subset 
U <Z N such that a\u = 0. This contradicts the statement that S{a) = N. Hence u is 
tangential to l. ■ 

An initial hypersurface associated with f is a submanifold a: S M such that 
each integral curve of v intersects S exactly once, and v contains no closed curves. 
For the following we further demand that each integral curve of v has domain M. If 
T is the parameter along an integral curve 7a;: M — ^ M, with 7a; (0) = x G S, then 
'^liAr) = Ix-kidrlr) = 7x(t). Thcsc two couditious imply that there is a diffeomorphism 
c: M X S ^ M with c(r, x) = j^^t). 

For M = £ and v = W, the Liouville vector field (|2]), then since W is horizontal 
the requirement that the domain of is M is equivalent to demanding that the proper 
time for each curve where = jx is from — oo to +oo. 

A box s: S M is any compact submanifold (usually with boundary) of an initial 
hypersurface and of the same dimension as the hypersurface. That is v is transverse to 
s and dimS* = dimM — 1. Given a box s: 5* ^ M we can define the expectation of 
as 

bW] = (^^^)[l] (61) 

which is valid since 1 G FAS' has compact support on S. 

The future of a box or an initial hypersurface is the set of points which lie on the 
integral curves of v in the direction of v away from 5* or S. 

The total charge on an initial hypersurface a: S "-^ M associated with \E' is given 

by 

Qtotai = (62) 

However in order to define charge we require that given a series of boxes Si C Sj+i 
with IJ S'j = E then limj^oof^i (^1^)] is well defined and independent of the choice of 

{^1,52,...}. 

Theorem 22 Given a solution a^a G FsAM to the transport equations where a: N 
M is an embedding. Let si. Si ^ M and S2'- S2 ^ M be two boxes with the property 
that any integral curve of v passing through one also passes through the other. Then 

[slia^a]] = [4(a^")] (63) 

proof: Assume first that 5*2 hes to the future of Si. Let r.U M be the closed 
submanifold with dimf/ = dimM, given by the union of all the integral curves of v 
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between 5*1 and 5*2. The pullback manifold with respect to l and a is given by a~^{U) 



and (49) becomes 


(defining 


the embedding 








(f/)V^iv 








a 




By definition s'i{a^a) = {a.i)^{si*a) so 

[s^(a^a)] = s^(a5a)[l] = (ai)5(si*a)[l] = D{s*a)[l] = j 



Si a 



a-^is,) 



Now a G TAN so i*a G rA(a ^U) is a regular form. Therefore we can perform 



the standard analysis. From the transport equations then (58) implies da = and (59) 
imphes iu^ = where u G TTN is the unique vector field satisfying aj,('u) = f |a(Af)• 
The boundary of U is given by 

d{U) = S2-Si + V 

where the — sign refers to the orientation of Si and u is tangential to V. Since a is an 
embedding the boundary of a^^{U) is given by 

d{a-\U)) = a~\S2) - a-\Si) + a-\V) 

where u is tangential to a^^iV). Therefore /g-ij-y) a = 0. Thus 

da = I a — I a + I a = [s{{a^a)] — [s2(a5a)] 







a-^u 



a — I a + I a 

a-i(S2) Ja-^Si) Ja-^v) 



For two general boxes Si and 5*2 simply take another box Ss lying to the future of 
both El and S2. ■ 

Corollary 23 Conservation of charge. Given two initial hypersurfaces a: S > M and 
a:fl ^ M and given that the charge (62) is well defined with respect to S then the 
charge is well defined with respect to S and 

[a^m = (64) 
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proof: Any box s: S corresponds to a box s: S > S. Therefore from theorem 

[s'^(\l/)] = [s'^(\E')]. Taking the limit of boxes shows that [cr'^(\E')] is well defined and that 



(64) holds. 



Given a submanifold distribution on an initial hypersurface we can generate a 
unique solution to the transport equations. 
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Theorem 24 Let a:J]'—>-Mbean initial hypersurface of v E TTM and /et c: M x S — >■ 
M he the corresponding diffeomorphism. Let b^a G TsAS be a submanifold distribution, 
with h:N ^T. andae rA'^'^^^A^. Let 

^ = a^(7r^a) G rsA'^""*^-'M (65) 

where 

a:Rx N ^ M ; a(r, y) = c(r, b{y)) (66) 
and 712'.^ X N N is the natural projection. Then \E' satisfies the transport equations 



(57) and the initial conditions a''^ = b^a. 

Furthermore \1/ is unique in that given any $ G FsAM which also satisfies the 
transport equations, and 0"''$ = b^a, then \E' = $. 

proof: To establish that \1/ is indeed a submanifold distribution, we must show the four 



conditions (41 44). The following commutes 
M X 



Idxb 




X S 



Thus a is an embedding. Since M x b{N) C M x S is a closed submanifold. Then 
a(A^ X M) c M is closed. Clearly S{7r*a) = M x A^. Also since Idxb:MxA^^Mx b{N) 
is a diffeomorphism so is a. 

To establish that d'^ = it is necessary to realise that M x A^ may have a boundary. 
Let l: B ^ N he the boundary of A^. Then the boundary of M x iV is given by 
hR X B ^ X N where t = Id x t, i.e. t(r, y) = (r, L.{y)). 

d"^ = da^{7T2<y) = da^Dijil^a) = a^dD{TT2<y) = a^{dTT2a) — a^L^{L*7T2a) 

Now (iTTgO; = 7T2da = since deg(a) = dimA^. 

Let Tf2'-R X B ^ B he the second projection so that 7T2 o t = l o 7:2. Then 





{tt2 o tya = {lo Tf2ya = 7r2(''*«) 



since dega = dimA^ > dim 5. Thus d'^ = 0. 
We now establish that i^^^ = 0. 



iya^[TC2a) 



i^a^D{-K*2a) = a^id^D{-Kla) = a^{iQ^7r2a) = a^iniKi.dr'^) = 



since a^,{dr) = v and 712*9,- = 0. 

We now establish that \1/ satisfies the initial conditions. Let lo'- N ^ M x A^, 
/o(?/) = (0,?/) then the following commutes and is the pullback manifold 



lo 



X N 



M 
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We now establish that \E' is unique. We first show that = S{^). Let s: S M 

be a box so that x G S" and fl S = 0. Now there hes a box s: S* ^ S which hes to the 



future or past of 5* hence from lemmas 22 and 18 we have 



Thus given x G then for any box s: S "-^ M with x G s{S) then = 

7^ 0. Thus X G «S(\I/) and visa versa. 
Since «S(\E') = S{^) we may consider up to diffeomorphism that $ = a^f3 for some 



/9 G rA(]R X N). From theorem 21, both vTga and /5 satisfy the transport equations 



with respect to 9,- on M x A^. Furthermore b^{lQTT2a) = a''^ = o"'^$ = b^{lQ(3) hence 

(M X A^). Hence 



IqP = Iq712(^, i.e. they agree on the initial hypersurface Iq: N 
j3 = TT^a and so $ = \E'. ■ 
There are two alternative formulations of one simply using the puUback, a second 



using an integral (65) which can easily be generalised for general distributions. 

M, b^a G FgAS, b: N 



Corollary 25 Let a-.T. ^ M , v e TTM, c: M x S - 
a G Vh!^''^^N and ^ G FsA'^^'^^-iM he as m theorem 

^ = (7r2 o c"i)^(b,a) 

where : M x S ^ S zs the second projection. 



Then 



(67) 



proof: Since c: M x S ^ M is a diffeomorphism then [7^2 o c ^): M — S. Thus to define 
the puUback [1x2 o c~^Y we have the following commutative diagram. 



X N 



M 



■K2 



N 



TT20C 



Clearly dim(]R x A^) + dim S = (dim A^ + 1) + (dim M - 1) = dim A^ + dim M so 712 o c^^ 



is transverse to b, and (67) follows from ^65^ and (51) 



Corollary 26 Let o-.H ^ M , v e TTM, c: M x S M, b^a E FgAS, b: A^ 
a G TA'^^'^^N and ^ G FsA'^''^^^-iM be as m theoremWA Then 



for all G FqA^M where $ = b,;Q; and Ct : S M, (riv) = cij^y). 



(68) 



proof: The pushforward of the vector 8^ \ i^r,y) G Ti^r^y) x N) under M x A^ ^-^^ M x S M 

is given by (Id X b)^{dr\{r,y)) = dr\(r,b{v)) cAlirMy)) = ^|c(T,b(y))- 

Let Ir : N ^Rx N, Ir{y) = (r,?/) so that (co (Id x b) o/^)(z/) = c( (Id x b)(r, y)) = 
c(i", b{y)) = Crib{y)). Thus 



a^{7!-2a)[(f>] = D{7i2a)[a* 



a (f) A TTjO; 



(Id X b)*(c» A 7r*a 



xN 
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dr I /*2ar ((Id x b)*(c» A vr^a) = f dr [ I* ({Id x b)^za,(c» A 7r*a) 



dr / I* ({Id X b)*(c%0) A 7i*a) = dr {co{ldx b)o Q*{iv(f)) A a 
I Jn ^ ' JN 



dr b*(CM))A«= / dTb,{a)[C{tv(p)] = / dr^Cr{^^ 



3.1. Distributional solutions to the Liouville equations as a source for Maxwell's 
equations. 

The following shows that the source for Maxwell's equations is closed. 

Theorem 27 The source for Maxwell's equations is closed i.e. dj = 0. 

proof: If G is bounded with respect to vr then lemma [6] implies (i(7r<^0) = 7r^{dQ). If B is 
regular and unbounded but n^Q is defined then in [3] it is also shown that don^ = n^od. 
Thus 

dj = d{7r,e) = 7r,{de) = 

■ 

We have established that a submanifold distributional solution to the transport 
equations can be constructed from a distribution on an initial hypersurface. From now 



on we consider distributional solutions for the transport equations (57) for distributions 
\E' G FsAii^ where n: E —>■ M is a fibre bundle. If this solution is to be the source for an 
electromagnetic field then we need to be able to take its pushforward with respect to 
the projection map vr. However this projection map is not proper. Therefore we must 
establish under what condition we can guarantee that \E' is bounded with respect to vr. 

Lemma 28 Let a: N ^ E and a^{a) G T^AE be bounded with respect to n: E ^ M 
then D{q) is bounded with respect to {tt o a): N M and 

TT^{a^a) = (tt o a)^a (69) 

proof: Let U G M he compact. Since a^{a) is bounded with respect to vr, it follows that 
Ti^^U n S{a^a) is compact. As a is proper a~^(7r^^f/ fl a(A^))) C is compact. Now 

a-\TT-^U n a(A^)) = {a'^TT-^U) n {a-^a{N)) = (vr o a)-^U f] N = {tt o a)-^U n S{Da) 

Thus D{a) is bounded with respect to (vro a). Given G FqAM and h G VqKE satisfies 



(34) i.e S{l — h) n«S(a*0) n«S(a^a) = then we wish to show that a*/i G FqAA^ satisfies 
(34) for \Ef = Da and replaced by 7r*0: 

5(1 - a*h) n <S(a*7r» n S{Da) 

= 5(1 - a*h) n 5(aV» nN = S{a*{l - h)) n 5(aV*0) 

C a-^(<S(l - h)) n a-^(<S(7r») = a-\S{l - h) n 5(7r») 

= a-\S{l -h)n S{7T*(f)) n a(A^)) = a-\S{l - h) f] 5(a» n S{a^a)) = 
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Hence we have 

7r<j(a^a)[0] = a^a[h7r*(j)] = Da[a*{hn*(f))] = Da[a*{h) {a*7r* 
= Da[a*{h) (vr o a)''(0)] = (tt o a)^a[(j)] 



giving (69). 



Recall that a vector field W G TTE with tt: E ^ M and E C TM is horizontal if 
7rv,(Vr|M) = u for all u & E. The following lemma relates the integral curves of horizontal 
vector fields with curves on the base space M. 

Lemma 29 Let n: E M with E C TM and W G TTE be a horizontal vector field 
and let 7: M ^ be an integral curve of W then 7 = C where C = vr o 7: M ^ M. 

proof: Let r G M. Recall that dr G T^-IR is the natural vector field and drlr G FT^M 
is a vector at the point r G M. The lifts of the curves 7 = 7^(9,-): M TE and 
C = a(5,):M ^ TM. 

C{t) = C,{dr\r) = (vr O ^)4dr\r) = TT^^Mlr)) = = MW^(r)) = 7^ 

■ 

We can ask the question: 

If given an initial hypersurface a:Il "—>■ E, with respect to a 
horizontal vector field, such that a'''^ is bounded with respect (70) 
to TT o o" then is bounded with respect to vr? 

It turns out that even for simple cases this is not the case. For example consider 
the Vlasov operator corresponding to a simple force free drift in one dimensional 
Galilean spacetime M = coordinates {t, x) with E = M.^ coordinated by {t, x, x) 
and 7i(t,x,x) = (t,x). 

A solution to the transport equations for this vector field is given by = a^a G TgA'^E 
where a: —>■ E, a(r, y) = (r, ry, y) and a = a{y)dy G FA-'^M^. 
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Now to see that \1/ is a solution to the transport equations we use theorem 21 
Clearly da = and 



f 



d 



d 



d 



Hence it is a solution. 

With initial hypersurface a:T, ^ E, cr(xo,Xo) = (^0^X0,^0) an initial distribution 
is cr''\l/ = b^a where b:M — S, b{y) = (toy,y), so that a{b{y)) = (to, toy, y). Given 
compact f/ C M we have 

Sia'"^) n n-^U = S{b^a) n n-^U = {{to,toy,y)\y eR} H {{t,x,x)\{t,x) E U} 

= {ito,toy,y)\ito,toy) eU} 

is compact as long as to 7^ 0. Thus is bounded with respect to vr o a, see figure [2] 
However consider (0, 0) E M then 

^"'{(0,0)} = {(0,0,a;)|x eR} = {{0,0,y)\y eR} = {a{0,y)\y eR} C S{^) 

Thus \E' is not bounded with respect to tt. 

It turns out that for the relativistic Maxwell-Vlasov equations, as long as the initial 
hypersurface S C is such that vr(S) C M is a Cauchy surface, then we can guarantee 
that if cr'"\E' is bounded with respect to tt o a then \1/ is bounded with respect to tt. 

Theorem 30 Let tt: S ^ Ai be the upper unit hyperboloid bundle over a globally 
hyperbolic spacetime. Let W E TTS be a horizontal vector field, i.e. such that 
T^*(W\u) = u. Let cx-.H "-^ S be an initial hypersurface such that 7r(cr(S)) G A4 is a 
Cauchy surface. Let \E' E FsA^^ be a solution to the transport equations with respect to 
W. Let cr'"(\E') be bounded with respect to tt o a then is bounded with respect to tt. 

proof: Let S = 7r(o"S), a: S ^ 7W and tt = tt o a : T, T,. Given a compact subset 
Y G A4, since S is a Cauchy surface, it is shown in |7] that for each y E Y the light 
cone of y intersecting with S is compact. Define X C S as 



X 



X G S 



there exists a timelike or lightlike curve C passing though x and Y 



then X is compact. Furthermore since X and Y are compact the proper time it takes 
to go from X to y has a maximum and minimum. 



T 



T E 



Then T C 



there exists a timelike normalised curve C : M 
such that C(0) E a{X) and C(r) E Y 

is compact. We have the following commutative diagram 

lo 



M 



X N 



M 
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Figure 3. On the lightlikc initial hypcrsurfacc, all initial velocities point towards the 
point (0,1) 



Since (t''\1/ = h^a is bounded with respect to tt, and X C S is compact then 
7T~'^{X) n S{b^a) C S is compact. Now 7i~\X) n S{b^a) = 7r"^(X) n b{N) and 
since b: N —>■ J] is proper 

b-\7T~\X) n b{N)) = b"^TT-\X) nN = b~'^TT-\X) CN is compact 

We need to show that 7i^^{Y) fl «S(\I/) is compact. Now from set theory and the 
fact that a is injective and <S(^) = a(A^) then 7r-^{Y) n «S(^) = a{a-^iT-^{Y)). Now 
since a is continuous, it is sufficient to show that a^^7i^^(Y) C M x is compact. 

Let {ro,y) G a~^iT'^(Y), so that 7r(a(ro,?/)) G Y. Let 7: M £^ be given by 



7(r) = a(r, y), then from (66) 7 is an integral curve of W. From lemma 29, 7 = C 



where C = vr o 7. Since C:]R — g{C,C) = —1 so C is a timelike normalised curve. 
Now 

C(0) = 7r(7(0)) = 7r(a(0,y)) = 7r(/o(2/)) = 7r(a(b(y))) = a(7r(b(y))) 

thus n{b{y)) G o-(S). Also C(ro) = 7r(7(ro)) = 7r(a(ro,|/)) G Y. From the definitions of 
X and T, C(0) G a{X) and tq G T. Since C(0) = <j(7r(b(?/))) G a{X) then 7r(b(?/)) G X 
so ?/ G b~^7r~^(X). Hence {TQ,y) G T x b~^7r^^(X). Thus we have shown that 

a-^-K-^Y) C T X b~^7r-^(X) 

Being the product of two compact sects T x b"^7r~^(X) is compact and since a~^7r~^(y) 
is closed it is therefore compact. Hence result. ■ 

The condition that 7r((T(S)) = S C is a Cauchy surface is required. Here we give 
an example where A4 is globally hyperbolic, ct'^(\E') is bounded with respect to n: S ^ M. 
but is not bounded with respect to tt. 

Consider Ai is two dimensional Minkowski spacetime, with coordinates (x°, x^) and 
£ is the upper unit hyperboloid coordinates {x^,x^,y). Let W be the Liouville vector 
(§ for a zero electric force: 

d d 



W = ^l + y^— + y 



dx^ 
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Let 

so that 71 o a-.H ^ Ai, 71 o a^z, p) = {z, z) is a lightlike curve so must be transverse 
to W . Also every integral curve of W must intersect cr(S) so S is an initial hypersurface 
of W . Clearly 7r(crS) is not a Cauchy surface. Let = M x G M, -S < 0} and 



1 — z ^ ^ —z ^ —z 



a: N ^ £ , a(r, i)=(x = =r + z,x = =r + z,y 



VI -2z ' VI -2z VI -2z 

so that 

/ d \ ( \ — z —z 



Hence a is tangential to W . Let \E' = a^l. Now (t'^(\I') = b^l where 

b: {5 G M, 5 < 0} ^ S , h{z) = {z, z, -z{l - 22)-^^^) 

This map is proper and cr'^(\E') is bounded with respect to vr. However this example has 
been contrived so that every integral curve passes through = (1,0) G Ai. See 

figure [3j 

a{Vl - 2z, z) = (1, 0, -z{l - 2z)~i/2) 

So 

7r-i{(l,0)}n5(^) = {(l,0,i/)|i/>0} 
is not compact so "if is not bounded with respect to tt. 

4. Known solutions to the Maxwell-Vlasov equations. 

4.I. Worldline of a point charge and the Klimontovich distribution 



For the worldline of a point charge we can show that the Liouville equation ( T5p6 ) 



implies that the point charge undergoes the Lorentz force equation (|3j) and that the 
electromagnetic source due to a point charge is given by ^7 = G FsA^A^ where 
= C^(l). However due to the divergence of the electromagnetic field F due to 
point sources we cannot demand that there exist solutions to the complete Maxwell- 



Vlasov equation (13 16). It is therefore necessary to perform a regularisation of F 
before substituting into For a single point charge the usual regularisation is the 
Dirac regularisation which lads to the Lorentz-Dirac equation. 

Let C:R ^ M he a regular curve with g{C,C) = —1 and let C:R ^ £ he the 
corresponding lift. Let 

e G FdA^^ ; Q = C^l = C^ i.e. e[0] = [ C*(0) for all G FqA^^ (72) 

Jr 

Clearly dQ = 0. 



Lemma 31 Equations (13) and (72) give the point source. 

J = C^ = (73) 
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proof: C is bounded with respect to vr thus from lemma 28 we have 
J = 7r,(0) = 7r,(C7^1) = (tt o C)^l = C^l = 



Lemma 32 Equations (16) and (72) imply that the point charge undergoes the 
Lorentz force equation: 



proof: From (16) and (72) we have iv/C = hence from lemma 



15 



(74) 

we have C is 
W\ 



tangential to W hence C is an integral curve of W with a scaling. I.e. w i^^-^-) 
k{t) C{t). By hitting both sides with vr^ we see 



hence k = 1. Thus W\ 



(74) 



C(r) 



C(r). So C is an integral curve of W and thus (3) gives 



As stated in the introduction the Klimontovich distribution is simply a finite 
number of point worldlines. Thus we can set 

N 
k=l 



Which gives the source for Maxwell's equation as 

N 



k=l 



Again, this does not give rise to a continuous electromagnetic field so some method 
of regularisation is necessary. For the Klimontovich distribution, for example, we can 
ignore the contribution of Cj to J when calculating the motion of the worldline Cj. That 
is we solve d{DFj) = and d{D(kFj)) = Ylik^tj^k- then construct the Liouville 
vector field Wj = W{Fj) using (b| and solve the transport equations for Cj using Wj. 



4-2. Cold charged fluid model 

Here we see that the cold charged fluid model[8] is simply an example of a distributional 
solution to the Maxwell-Vlasov equations. Let v G TTAi with g{v,v) = —1 and 
p G TA^Ai. We may consider f: £^ as a closed embedding. Let 

v,{p^l^ (75) 
Q\(i)]= I ^\(vL) Ap\x*v\x for all ^gFoA^^ 



6 G FnA^^ 



i.e. 



e 



where v G FA^A^ is the metric dual of v satisfying v{u) = g{v,u) for all u G TTAi. 



Lemma 33 Equations [T^ and (75) give the source for the cold charged fluid: 
J = D{p^v) 



(76) 
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proof: is bounded with respect to tt thus from lemma 28 we have 
J' = n^Q = 7r^(f^(p * v)) = (tt o v)^{p •kv) = D{p -k v) 



Also since dj' = we have the continuity equation 

d{p i^v) =0 



Lemma 34 Equations [^,(16) and (75) imply the Lorentz force equation 



(77) 



(7J 



proof: From (16) and (75) we have iw{v^{p -kv)) = 0. Thus from lemma 15 there exists 
u G TTM. such that Vi,{u\x) = W\y\^. Thus u\x = ir^Vi.^ulx) = vrj,H^|„|^ = v\x hence 
u = v. Hence f*(f = W\^\^. 

Given an integral curve C: M ^ oi v then f |c(t) = C{t), i.e. v o C = C. Thus 

= M^lciro)) = W"l.|c(.o) = ^lc(ro) 

Hence C is an integral curve of W and hence from (3) VqC = i^F. Since this is true 



for all integral curves we have (78) 



4-.3. Multicurrent model 

Here we see that the multicurrent model is simply an example of a distributional solution 
to the Maxwell-Vlasov equations. The conclusions ( [8l|82[84[ ) are in [9l|10], with slightly 
modified notation. An application to wake breaking in plasmas is given in [TT]. 

Let i3 = M X ^ be a body-time manifold with points (t, G B and a measure 
rfr A /C G TA^B where /C = 7r^(/C) G TA^B and IC G FA^^.'het C:B M with 
C = C^dr) G S. Let 

e G FdA^^ ; e = Q(/C) (79) 



I.e. 



&[<!>]=[ 0b(ry)A/C for all ^gFoA^^ 

J{T,y)eB 



Lemma 35 Equations (13) and (79) give the source J due to a multicurrent fluid. 
These may be written in a variety of ways. Using distributional push forwards: 

J = Q(/C) (80) 

Furthermore if J is regular in the sense that it is piecewise continuous, then we can 
write J = D{J) where J G T^^A'^M. is piecewise continuous. 
In terms of the inverse pull back on generic open sets 

J= sign(det(q,].))q^i'^(/C) (81) 

1=1 
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where t/^ G Ai is an open set where C is generic and therefore the number of preimages 
f/j^, i = 1, . . . , N{U^) is constant, and the maps C[j]: f/j^ — > are invertible. In terms 
of integrals over spacelike hypersurfaces of Ai. 



J 



/C for all bounded spacelike hypersurfaces S G Ai 



proof: Since C{T,y) = C^{dr\(r,y)), 

T^{C{T,y)) = 'KC^{dr\(r,y)) = C{T,y) 

i.e. IT o C = C . Hence from lemma [isl 

J = 7r,(Q/C) = (tt o C),{}C) = Q(/C) 



(82) 



giving (80). 

Let G M he a generic subset and G TqA^U^ then since C([/|) = ±U^ 
depending on the sign(det(C[i]^)) we have 



A J = D{JM = J[ct>] = C,{IC) 



M 



Oa/C 



N{U-^) 

E 

i=l 



Cr1,(/'A/C 



2=1 



Si: 



ign(det(q,].))q^^^(q]0A/C) 



sign(det(q,],))</)Aq^^*(/C) 



since this is true for all (p G FoA^f/^ then we have (81). 

Putting s:S'-^Ai and C: i3 ^ into (|48]) we see that {(y, x) G i3 x S" | C{y) 
s{x)} = C~^{S) thus (49) becomes 



c 

Y 

— 



c 
M 



(83) 



where for y G C ^{S) we have C{y) = C{y) and s{y) = y. Hence from (80) and since S 
is compact so 1 G FoAS* we have 



s*J = D{s*J)[l] = s'{DJ)[l] = s'J[l] = s'{C,IC)[l] = C,{s*IC)[l] 



I.e. (82). 



Lemma 36 Equation [^,(16) and (79) imply 



15 



(84) 

there exists 



proof: From (16) and (79) we have iwC^{IC) = 0. Thus from lemma 

U G TTB such that Ci,{u\(^r,y)) = ^^^(rj/)- 

Hence u = dr- Thus C^{dT-\(T,y)) = W\c(Ty)- Keeping y as a constant then C{T,y) is 



an 



integral curve of W and hence form ^ we have (84). 
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4.4- The water bag model 

Let C be a bounded 3 dimensional manifold with boundary l: B ^ N . Let 
a: X ^ £^ be an embedding such that o a = 7^2'- N ^ M. ^ M. is the second 
projection. This induces the map v = {a o t): B x Ai ^ S so tt o v = tt2- Let 
a = a*{iw^) G rA'^(iV x M). We call this the water bag model 0. 

The source for Maxwell's equations for the water bag model is given by 

J' = vr^a^a 

Comparing with (ml), this is a regular distribution J' = D{-kJ) where j E TA^A4 



^ J7r-i{x)na(Af) yo ' 
Let V = ao l: B X J\A ^ £ and v{^, x) = v^{x) so for each ^ E B, E VS. There is 



no requirement that satisfies the cold charged fluid equations (78). However we can 
find a reparametrisation = W(^(^^^x){x) = v^{x) such that does satisfy (78). 

To see this, let S G be a Cauchy hypersurface so that tt^^S C £^ is an initial 
hypersurface with respect to W. Thus we have the submanifold a: B x H ^ S. Given 
any {(,y) E B xT, let Co-{c,y)-^ — > £^ be the integral curve of W passing though a{(, y), 



which therefore satisfies (74). Given ^ E B and x E Ai, since v^{x) E v{B x A4) then 
v^{x) lies on some integral curve Ca-{(,y) and hence there exists unique r E M., ( E B 
and y E such that v^{x) = C„((^^y){T). Let C{^^^) be the corresponding ( so that 



v^{x) = C'^(c(c,x),j/)(t")- Now set = v^{x) so obeys ( |78| ). 

5. Conclusion and Discussion 

We have presented the distributional Maxwell-Vlasov equations and analysed the 
consequences of pushforward distributional solutions. We have paid particular attention 
to the question of when a pushforward distributional solution can be the source for the 
Maxwell equations. 

There are a number of directions that this research may continue: 
It is possible to extend the idea of a submanifold distribution to include distribution 
with support on submanifolds with different dimension or intersecting components. One 
can ask how to extend the notion of a puUback to these distributions. From equation 



(68) we can see it is possible to take any distribution on an initial hypersurface and 
create a distributional solution to the transport equations. 

Submanifold solutions with dimension of 5 or 6 may, for example, correspond to 
Vlasov one particle probability functions which at every point in spacetime have a 
range of possible values of one component of velocity, but fixed values for the other 
components. 

As stated in the introduction, it may be possible to solve the Maxwell-Vlasov 
equation numerically, for distributional solutions of arbitrary dimension. For low 
dimension it will be necessary to perform some form of regularisation of the 
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electromagnetic field F, however this should be easier for the solutions of dimension 
2 or 3. These numerical solutions can either be interpreted in their own right as, for 
example, the dynamics of strings or disks, or as approximations to regular solutions. 

The stress energy tensor for a distributional source on the upper unit hyperboloid 
is known. This enables one to couple the Maxwell-Vlasov equations with Einstein's 
equations to give the Einstein-Maxwell- Vlasov system [12]. Since the equations in this 
article are already covariant no further modification would be required. 
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